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The goal of this study is to investigate the scaling behaviour of our 2 HEX action. For this purpose, 
we compute the Nf = 3 spectrum and compare the results to our 6 EXP action. We find a large 
scaling window up to ~ 0.15fm along with small scaling corrections at the 2%-level and full 
compatibility with our previous study. As a second important observable to be tested for scaling, 
we chose the non-perturbatively renormalized quenched strange quark mass. Here we find a fairly 
flat scaling with a broad scaling range up to ~ 0.15fm and perfect agreement with the literature. 
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1. Definition of the 2 HEX action 

For the gauge part in the Nf = 3 simulations, we use a tree-level improved Symanzik gauge 
action 



j ReTr ( 1 " ^Plaq) + J ReTl " ( 1 " U « 
plaq rect 



(1.1) 



where c\ = —1/12 and cq= 1 — 8ci = 5/3 and the standard Wilson plaquette action in case of the 
quenched study. For the fermionic part, we used a tree-level improved clover-action 



Si = Sw ~ ^ Yl E (^°Mv^v [V] V) (*) > (1-2) 

X /i<V 



with csw = 1 and Sw the standard Wilson action, coupled to smeared links V. These smeared links 
are constructed by combining the HYP setup [1] with the analytic stout (EXP) recipe [2]: 

r&pW = U a (x)U^x + a)ut(x + fL) 

V@p(x) = e X p(^P TA {rJ i 1 : ) vp (x)^( x )})^ (x ) 
lf v (x) = Y ^ v (x)V'(^(x + a)vS t v (x + A) 

V®(x) = e X p(^P TA {r( 2) v (x)^(x)})^(x) 

±v^u 

V M (x) = e X p(^P TA {rt 5) (x)[/;(x)})[/ M (x). (1.3) 

This smearing was introduced in [3]. Here we chose the weights (oti, CXq., CC3) = (0.95,0.76,0.38). 
It is straightforward to show that this smearing is analytic with respect to the thin links and 
hence can be used within an HMC for dynamical simulations. We applied the recipe (1.3) twice 
to define our 2 HEX action. Generically, such smearing improves the chirality of the underlying 
Dirac operator and drives renormalization constants closer to their tree-level values [3]. Hence we 
expect a scaling behaviour which is close to that of a non-perturbatively 0{a) improved action, 
although our action is formally G(ad) improved. 



2. Scaling of Nf = 3 hadron masses 

As a first test for this action, we computed the Nf = 3 hadron spectrum to be able to compare 
the results to our former 6 EXP action [4]. We took four values for j8 so that the cutoff varied 
between a w 0.06 fm and 0.2 fm. At each beta we simulated at least four different masses to be able 
to perform a safe interpolation to our reference point M n /M p = [2(M£ hys ) 2 - (M£ hys ) 2 ] 1/2 /Mf ys 
0.67. The hadron masses themselves were obtained by applying correlated cosh or sinh-fits to the 
correlators, where we avoided fitting excited state contributions by inspecting the effective mass 
plateau and chose the fit ranges accordingly. The PCAC mass was computed by fitting the plateau 
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of (doAoP) I (PP). Then we extrapolated the octet and decouplet masses to the continuum assuming 
either G (ad) or G (a 2 ) scaling, as displayed in Fig. 1 . We see flat scaling with a maximum of ~ 2% 
scaling corrections for the delta mass around a ~ 0.16fm. The scaling window is broad and extends 
at least up to this lattice spacing. The fit qualities favor the G(a 2 ) -scaling but in order to compute 
a reliable systematical error, the G(aa)-sca\ing assumption should also be taken into account. 
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Figure 1: Continuum extrapolation of the octet and decouplet masses assuming different scaling. 



3. Scaling of quenched quark masses 

Another important observable to study is the strange quark mass. The renormalization factors 
have been computed by using the non-perturbative RI method [5]. The renormalization conditions 
are defined as follows: 

Z {na,g(a))Z Q \na,g(a))r (pa)\ p 2 =tl 2 = 1, (3.1) 

where Yo{pa) = Tr{(S~ l GoS~ l )(pa) Po] is the renormalized vertex function belonging to operator 
O with the unamputated Greens function Go and a projector Pq, which projects to the quantum 
numbers of O. In order to remove an G (a) contact term we applied the trace subtraction procedure 
described in [6-8], therefore replacing all quark propagators S by S = S — TrpS/4. The quark 
mass renormalization constant is l/Zf 1 , but the renormalization condition (3.1) only gives ratios 
(Zo/Zq) ri . In order to avoid discrete derivatives as they would appear in the calculation of the 
wave-function renormalization Zq itself within the RI formalism, Zq was computed indirectly by 
using the VWI. To be more precise, consider the ratio 

^ {} (P(T/2)P(0)) ' K ' 

with P(x,t) = Y^ x W\ 75^2^) t) and V^{x,t) = ^2 X y/2 74 V2 The tree-level on-shell-improved 
vector renormalization constant can then be obtained by computing [9, 10] 

Z V (1 +am w ) = [C(t > T/2) - £fo - T/2)]-\ (3.3) 
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where m w = m — m cnt . Hence the wave-function renormalization can be obtained via Zq = 
(Zq/Zv^Zv- We found this method to be more precise than changing to the RI' scheme or using 
the conserved vector current within the RI formalism. 

Our data are reasonably consistent with perturbation theory for p > 3GeV. We found that the 
influence of the cutoff on the RI data is small for momenta up to p < n/(2a), i.e. the scale p « 
3GeV was not safely reachable at the coarsest lattice as can be seen from Fig. 2. Hence we 
computed the ratio 

*( M V0 = limJ^^, (3.4) 
a->0 Zg(H ,a) 

for /i' = 3.5 GeV and \i" = 2.1GeV only on the three finest lattices. The continuum limit was taken 
again assuming ff(aa) or &(a 2 ) scaling and is relatively flat (c.f. Fig. 3), as suggested by Fig. 2. 
Using these relations, we defined the quenched quark mass via 

m vw V) = (l-am w /2)m w /[R{^',^')Z s (^",a)]. (3.5) 

The continuum extrapolations are displayed in Fig. 4. As before, we considered both versions 
of scaling and see that both extrapolations are compatible and fairly flat, with at most 7% scaling 
corrections up to ~ 0.12fm. When considering the full dataset, the 0(a 2 ) extrapolation is slightly 
favored. Our combined result is 

[(m s + m ud )r ] m (2GcV) = 0.2609(39) (28), (3.6) 

where the first error is statistical and the second error systematical, not including the uncertainty due 
to quenching. This is result is compatible with [11] and in good agreement with other results from 
the literature. Assuming m s /m uc i « 27.5 as suggested by recent unquenched lattice calculations 
and ro = 0.49 this translates into 

mp(2GeV) = 101 .4( 1 .5) (1 . 1 ) . (3.7) 

Note, that for this value there is an unknown systematic error due to the above assumptions for 
tris/mud and ro and the systematic error analysis described below does not include these sources of 
uncertainties. 

The statistical errors have been computed by using 2000 moving43lock43ootstrap samples with 
atomic blocking because all configurations are well decorrelated. For handling the systematic 
eiTors, we followed the procedure described in [12]. Thus we performed the analysis using three 
different fit ranges to extract the masses from the correlators, two different scaling assumptions 
for the continuum limit and three scales pL = 3.0, 3.5 and 4.0GeV at which we matched non- 
perturbative results onto perturbation theory. This yields 3-2-3 = 18 different continuum limits 
which were all weighted by their quality of fits. The mean gives the overall estimate and the 
variance our systematic error. 
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Figure 2: Multiplicatively matched Z5. The colored dashed vertical bars denote the different cutoff p 
%l (2a) for the three coarsest lattice spacings. For j3 = 6.3, (jJ.ro) 2 ~ 179 is off the scale. 
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Figure 3: Continuum limit of Z s (3 .5 GeV, a) /Z s (2. 1 GeV, a) . 



5 



Scaling study for 2 HEX smeared fermions 



Thorsten Kurth 



0.34 
0.32 
0.3 

O 
i_ 

+ 0.28 

"O 

E 

0.26 
0.24 
0.22 




0.34 



0.32 



0.3 



+ 0.28 



0.26 



0.24 
0.22 



Garden et al., [Nucl.Phys.B 2000] 
aa-extrapolation 



0.02 



0.04 



0.06 



0.08 



aa/r r 



(a) ^"(aa)-scaling 




a -extrapolation 



0.02 



0.04 



0.06 



0.08 



(a/r o r 



0.1 



Garden et al., [Nucl.Phys.B 2000] ■ — b- 



0.1 



(b) ^(a 2 )-scaling 

Figure 4: Continuum extrapolation of the quenched (m s +m uc i) m at 2GeV in units of r$. 
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4. Summary 

We performed hadron and quark masses scaling studies for our new 2 HEX action. We found 
large scaling windows and mild deviations from the continuum limit in both cases. Moreover, 
this action significantly reduces the condition number of the Wilson-Dirac operator for small quark 
masses. We conclude that this action should be well suited for large unquenched phenomenological 
studies. 
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